Variational Theory of Flux Line Liquids 
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We formulate a variational (Hartree like) description of flux line liquids which improves on the 
theory we developed in an earlier paper [A.M. Ettouhami, Phys. Rev. B 65, 134504 (2002)]. 
We derive, in particular, how the massive term confining the fluctuations of flux lines varies with 
temperature and show that this term vanishes at high enough temperatures where the vortices 
behave as freely fluctuating elastic lines. 
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I. INTRODUCTION 

The study of the properties of flux line solids and liq- 
uids in high temperature superconductors (HTSC) has 
been, during the past few years, one of the most ac- 
tive areas of research in the field of superconductivity. 
One theoretical approach that has had a rather strong 
impact on our present understanding of flux liquids in 
HTSC is the boson mapping 1-5 , which is based on the 
observation 6 that there is a formal mapping between the 
partition function of a three-dimensional system of inter- 
acting flux lines, and the imaginary-time partition func- 
tion of quantum bosons in two-dimensions. One of the 
main results of this approach, which has in fact been 
reproduced using other methods 7 , is that the structure 
factor 



S(r,z) = (p(r,z)p(0,0)) 



(1.1) 



of a liquid of flux lines described by the trajectories 
Ri(z) = (vi{z),z) as they traverse the superconducting 

sample (here p(r,z) = Yli=i — r i( z )) i s the density 
operator in the r = [x, y) plane at heigh z, N being the 
total number of vortices in the sample), is such that the 
partial Fourier transform S(q, z) — (/3(q, z)p{— q, 0)) is 
given by 



S(q,*) = S(q,* = 0)e- 



<q)\*\/T 



(1.2) 



where e{q) is the excitation spectrum of the equivalent 
boson superfluid (whose precise form will be given be- 
low), and T is temperature (throghout this paper, we use 
units such that Boltzmann's constant ks — 1). Although 
the Boson mapping (and other similar approaches 7 which 
use the density as the basic dynamical variable) does not 
make any explicit predictions for the correlations of con- 
formational variables of flux lines, the exponential decay 
(1.2) of density correlations along the direction of the 
lines is generally interpreted as an indication that the 
flux line liquid is in a heavily entangled or "superfluid" 
state characterized by correlations of the form 



with a "diffusion" constant D of order (here n is the tilt 
modulus of the flux lines) : 



D 



T 



and by a mean projected area of the flux lines : 



DL 



u(z) = r(z) - (r(z)) 



(1.4) 



(1.5) 



which diverges with the superconducting sample thick- 
ness L. 

Since its inception by Nelson and Seung 2 , the boson 
mapping has had a widespread acceptance and a rather 
strong impact on our understanding of flux line liquids. 
It has been applied, in particular, to disorder-free vortex 
liquids 5 , as well as to vortex liquids in presence of point 3 
and correlated 8 disorder. In all these studies, the boson 
mapping gave seemingly reasonable results, except at one 
occurrence where it was noticed by Tauber and Nelson 4 
that the boson mapping gives nonsensical results in the 
presence of splayed columnar disorder 9 . 

In a recent manuscript 10 , henceforth referred to as 
(I), I have proposed a new approach to study three- 
dimensional flux line liquids in type II superconductors. 
This approach, which makes contact with the standard 
theory of classical fluids, is based on the separation of the 
flux lines' conformal variables into center of mass (cm.) 
and internal modes, and on the observation that the in- 
teractions between flux lines may lead, under certain con- 
ditions, to the confinement of the internal modes whose 
fluctuations are now bounded and no-longer diverge with 
the sample thickness L as in Eq. (1.5). While, as men- 
tioned above, the Boson mapping does not deal with the 
conformation variables of flux lines themselves, the the- 
ory developed in (I) predicted that the effective Hamil- 
tonian of the internal modes of interacting flux lines in a 
vortex liquid has the form (i here labels flux lines, and 
N is the total number of vortices in the sample) 



N 



i=l 



dz 



1 fdUiV 2 1 o ] 
2 K Vdz-) + 2^ Ui{z) \ (L6) 



<|r( 2 )~r(0)| 2 )=2,D|z| 



(1.3) where the mass p,o is given by 



1 



Mo = f|d 2 rV 2 V(r). 9o (r). (1.7) 

In the above equation, V is the interaction potential 
between vortices (see below, Eq. (2.1)) and go is the 
pair distribution function of the two-dimensional liquid 
formed by the centers of mass of the flux lines. Because 
in (I) we only used a simple Taylor expansion of the total 
Hamiltonian of flux lines in terms of the Uj's, the "mass" 
term /xo does not fully capture the effect of thermal fluc- 
tuations, and varies with temperature only through the 
(weak) variation of go- In particular, Eq. (1.7) shows 
that the value of the "mass" term fi remains almost 
unchanged as temperature T — > oo, which of course is 
unrealistic, since we expect the vortex liquid to behave 
as an ideal gas of freely fluctuating flux lines in this limit. 

In the present paper, we generalize the method pro- 
posed in (I), and construct a variational theory of inter- 
acting flux lines in the vortex liquid state. This varia- 
tional method, which in fact turns out to be nothing but 
a self-consistent Hartree approximation, captures the ef- 
fect of thermal fluctuations in a better way and leads to 
a "renormalized" mass term /x(T) which vanishes at high 
enough temperatures, as suggested by physical intuition. 

Our method of approach will be as follows. In section 
II we construct our variational approach and derive the 
expression of the renormalized fi(T). In section III, we 
compare the method developed in (I) and in the present 
paper with the approaches of refs. 1-5 ' 7 , and discuss in 
particular the interacting structure factor S*(r, z) derived 
using both methods. Section IV contains our conclusions. 



II. VARIATIONAL THEORY OF FLUX LINE 
LIQUIDS 

We thus consider a liquid of interacting flux lines 
in a d = (dj_ + l)-dimcnsional superconducting sam- 
ple of thickness L. We here will use the following 
Hamiltonian 2 ' 8 

*{j«(^)'+5E "('.<•>-'-«)} 

(2.1) 

where n is the tilt modulus per unit length of the flux 
lines, and V(r) is the interaction potential between flux 
line elements at equal height z. For a uniaxial HTSC, 
with both the average vortex direction and the unit vec- 
tor z aligned with the c axis, k ~ e 2 £o and V(r) = 
2soK (r)/\ a b, where the energy scale £o is given by 
£o = ( ( / ) o/47rA 2 b ) 2 , e = A a t/A c is the ratio of the Lon- 
don penetration depths in the (ab) planes and in the 
direction of the c axis, respectively; (f>o = hc/2e is the 
flux quantum, and Kq is a modified Bessel function 18 . In 
the developments that will follow, we will find it useful 



to work in Fourier space, and to write down the follow- 
ing decomposition of the flux line trajectories r^(z, t) into 
Rouse 19 modes 

oo 

r i (z)= rMnV q " z (2.2) 

n— — oo 

where q n = 2nn/L, and where the Fourier coefficients 
ri(q n ) are given by 

r l (qn) = y [ L dzr^e"^- (2.3) 
L Jo 

It will also be convenient to write Yi{z) as the sum 

Ti(z) =r 0i (t)+Ui(z,t) (2.4) 

where i*oi(i) = Vi(q n = 0, t) is the cm. position and Uj(z) 
the flux line displacement with respect to the center of 
mass position at height z. 

The statistical mechanics of the flux line system, Eq. 
(2.1), is described by the partition function 

N 

Z= HldnWe-f" 1 (2.5) 

J i=i 

/N . N 

l[dv 0i / J][dui(^)]e-^ (2.6) 
i=l i=l 

where (3 = 1/T is the inverse temperature, and where in 
going from the first to the second line we separated the 
path integral over every flux line trajectory into a path 
integral over all internal modes and an ordinary inte- 
gral over the cm. coordinates. The integration measure 
[duj(z)] is given by 

oo 

[dUi(z)] = J| du re (q n ) du lm (q n ) (2.7) 

71=1 

with u re (q n ) and u im (q n ) the real and imaginary parts 
of u(q n ), respectively. Now, if we were able to perform 
the functional integrations over the internal modes, such 
an integration would give us the expression 

N 

Z= \\dr {H c -f )H ^f K"k}] (2.8) 
■* i=i 

which looks like the partition function of a system of 
N ordinary classical particles with spatial coordinates 
{r i}, interacting through the Hamiltonian H e ff[{r oi }], 
and whose thermodynamics can be studied by standard 
methods of statistical mechanics. Unfortunately, since 
the internal modes {ui(z)} appear explicitely in the ar- 
gument of the interaction potential V in Eq. (2.1), it is 
not possible, in general, to perform the functional inte- 
grations over these variables in closed form. In (I), we 
used a simple Taylor expansion of the interaction poten- 
tial V(n(z) - Tj(z)) = V(r i - r j + u t (z) - Uj(z)) in 
terms of the u,'s to be able to perform the integrations in 



2 



Eq. (2.6) and make predictions for physical observables 
of our flux liquid. Here, and in order to better capture 
the effect of thermal fluctuations, we want to improve 
on the above Taylor expansion and use a variational ap- 
proach to derive the statistical mechanics of our model 
system. 

In general, variational approaches to classical statistics 
are based on the Jensen-Peirels inequality 21 



Z > Z„ 



Z x e 



-0(H-H 1 ) 1 



(2.9) 



where Hx is any Hamiltonian, and where the av- 
erage is performed with the statistical weight 
exp(-/3Hx)/Zx, with Z x = Tr(c~^ ffl ). While an arbi- 
trary choice of Hx will not necessarily lead to good re- 
sults, a judicious choice of the variational Hamiltonian 
Hx (that is simple enough so that thermal averages can 
be calculated and at the same time general enough to 
capture the physics of the full Hamiltonian) can lead to 
a variational free energy 



F v = -T\rxZ v 

= -TlixZx + (H- Hx)x 



(2.10) 



which is a very good approximation of the true free en- 
ergy of the sytem, and to an accurate description of the 
overall behaviour of the original model. 

We now have to choose a variational Hamiltonian to 
approximate the Hamiltonian H of equation (2.1). In 
what follows, we shall assume that we can write for our 
variational Hamiltonian Hx a decomposition of the form : 



Hu[{iii(z)}] describes the elasticity and mutual interac- 
tions of the internal modes of flux lines. For H u [{ui(z)}], 
we shall use the following Gaussian approximation (here 
and in the following, we sum over repeated indices unless 
otherwise indicated) : 

H - = ^2ZZ J dz J dz ' [G- 1 (z-z , )] ij u i (z)-u j {z') 

(2.12) 

Throughout the rest of the paper, it will be understood 
that integrals over the z and z' variables run from to L. 
The above expression can be rewritten in Fourier space 
in the form 

H ^ = \J2J2 [G-'MkuM ■ Uji-qJ (2.13) 

i,j n^O 

In Eqs. (2.12)-(2.13), [G _1 (ft0]y is an N x N matrix of 
variational parameters describing the interactions of the 
internal modes of flux lines 23 . In the boson language, the 
decomposition (2.11) of the total Hamiltonian into cm. 
and internal modes pieces can be seen as the the gener- 
alization of the Feynman-Kleinert variational method 20 
for single quantum particles in imaginary time to an as- 
sembly of interacting quantum particles. Both quantities 
Hq and [G^ 1 (q n )]ij will be determined variationally, and 
one can in fact show, following ref. 28, that the Gaussian 
approximation (2.13), with the optimal choice for the 



Hx[{r 0t }, {u,(z)}] - H [{r 0l }} + H u [{u t (z)}} (2.11) propagator [G 1 ( g „)] y to be detemined below, becomes 



In the above decomposition, H [{r i}} = 
H (r 01 , . . . , r jv) is an effective variational Hamiltonian 
describing the interactions of the cm. modes, while 



exact in the limit of a large number of components of the 
displacement field {ui(z)}. 

Using the Hamiltonian Hx above, Eqs.(2.11)-(2.13), to 
evaluate the variational free energy in Eq. (2.10), we find 



F v = -TlxxZx + ^ [ f[dv 0i $(a)c-^[^^^TL([Go 1 ( g „)] y - - [G" 1 ^)^) Gy(«„) + 

1 i=X i,j n/0 

+ \ L E/ nq)c^ r ™)e-^ -H [{r 0i }]\ 



(2.14) 



where 



i,j n^O 



Zx= f f[dr (H J [dvaiz)} e-? H ° cxp ( - | ]T ^(G" 1 ^))^-^^) • u 3 (-q n )) 

i=X 
N 

= I l[dr (H $(G)e- 0H ° 



(2.15) 



and and where G(q n ) denotes the Nx N matrix dj (q n )- The 

correlation function dj is on the other hand given by 



« 5 » - ft (aasknf ■ «" 6 » 



1 



n{ dc t(Tg(q n ))J ' C« = — (M^-u^) (2.17) 
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Variation of the above expression with respect to the cm. 
Hamiltonian Ho shows that F v is minimal for the partic- 

1 



ular choice 



i.j n^O 



i,3 



(2.18) 



The trial free energy F v of Eq. (2.14) thus reduces to pression (2.18) of H should be used for H . 

F v ThiZi (2-19) Further minimization of F v with respect to Gij(q n ) 

where now in the expression of Zi, Eq. (2.15), the ex- leads to the following results (Appendix A) 



[G-HQnKv = ml ^^ 2 ^(q).9o(q)e- 



9 2 <« 2 >/cil 



for i = j 



d ± (N- 



T) J ?V(q) 5o (q)e- 92 < u2 >/^ for i+. 



(2.20a) 
(2.20b) 



where the tilde indicates that the inverse propagator 
G^ 1 (q n ) has been averaged over all possible configura- 
tions of the cm. positions {r i} compatible with a liquid 
structure. 

The integral on the right hand side of Eq. (2.20b) being 
finite and independent of the total number of vortices N, 
we see that the off-diagonal elements [G~ 1 (q n )] i ^ vanish 
in the thermodynamic N — > oo limit, in which case the 
internal modes elastic propagator is diagonal and given 
by: 



(2.21) 



with Sij the Kronecker delta symbol and where p is given 
by: 

M = -{^j/ V ^M)^ q2{u2)/dl - (2-22) 

In the above equation, the mean square width of a flux 
line 



(u 2 ) =d 1 _Tj2G(q n 
dj_T 



(2.23) 



depends on fi. Thus we see that Eq. (2.22) is in fact 
a self-consistent equation for the "mass" coefficient /i. 
From this last equation, we also see that /i depends on 
the density p of flux lines, not only through the pref- 
actor in front of the integral, but also through the pair 
distribution function <7o( r )- Through this last quantity, 
p also depends on temperature and on higher, nontriv- 
ial correlations between the positions of flux lines in the 
vortex liquid. The strongest dependence of p on tempera- 
ture, however, occurs through the "Debye- Waller" factor 
e -« ( u )/ d ± which was missing in our more elementary 
Taylor expansion in (I) , and which describes the effect of 
internal fluctuations of flux lines. 



From Eq. (2.22), it is not possible to derive a gen- 
eral expression of the coefficient p without making as- 
sumptions about the analytic form of the pair correla- 
tion function <7o( r )- Here, we shall place ourselves in the 
case of a dilute flux liquid and use the analytical ansatz 
we proposed in (I) for the pair correlation function go(r), 
namely : 



ffo (r) = 1 - i]cxp(-ar 2 /a 2 ) 



(2.24) 



where a is a constant of order unity, and < r\ < 1. The 
numerical constant rj quantifies the correlations between 
cm. positions of flux lines. It is close to unity when flux 
lines are strongly anti-correlated due to the repulsive in- 
teractions between their surrounding supercurrents, and 
close to zero in situations where there is considerable cut- 
ting and crossing of flux lines. Using the above ansatz for 
go(r), we obtain the following expression for the "mass" 
p of the internal modes as a function of T: 



KT) = Mo 




(2.25) 



where /xo = ^{T = 0) = 2rjnpeo/d±. Fig.l shows a plot of 
p(T)/po as a function of the variable v = aT / 4a 2 ^/npo . 

From Eq. (2.25), it is not difficult to see that the above 
expression has the limiting value n(T) — > it when T <C 
a 2 ^Kp~o 7 while at higher temperatures T 3> a 2 y/Kp , 
fi(T) ~ /io(4ct 2 y/ K/io / aT) and goes to zero like \/T. We 
thus see that (u 2 ) in Eq. (2.23) will vary with tempera- 
ture like T 3 / 2 . More generally, any nontrivial form of the 
pair distribution function go(r) in Eqs. (2.20a)-(2.20b) 
will result in a nontrivial temperature dependence of 
/i(T) and hence of the mean projected area of a flux line 
(u 2 ) . This temperature variation of (it 2 ) can be measured 
in numerical simulations and is actually a good way to 
test the predictions of the present paper. Note also that 
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FIG. 1. Plot of n(T)/no versus v = aT / '4a 2 ^/np^ . 

as T becomes very large, p{T) vanishes and we recover a 
liquid of freely fluctuating flux lines with (it 2 ) diverging 
with the sample thickness L. 



ILL THE INTERACTING STRUCTURE FACTOR 
OF FLUX LINE LIQUIDS 

We now turn our attention to a comparison and a dis- 
cussion of the results obtained by the approach developed 
in (I) and the present paper, and the Boson-mapping type 
of approach. Since this last type of approach only makes 
predictions for correlations of coarse-grained variables, 
such as the density, and makes no predictions for micro- 
scopic conformation variables of the flux lines, we here 
shall focus on comparing the interacting structure factor 
S(r, z) obtained through both methods. 

Let us first start by deriving the structure factor of a 
liquid of flux lines from our variational approach. This 
can be done in a straightforward way, with the result 
(Appendix B) : 



5(q,z) =p 2 .9o(q) 



2 (u 2 )/d A 



pe 



-q 2 <P(z) 



(3.1) 



Let us insist that this result is actually very general, and 
only requires that an approximate decomposition of the 
form (2.11) can be written for the Hamiltonian of the flux 
lines in terms of the conformational degrees of freedom 
of flux lines, and that the internal modes of flux lines 
obey a Gaussian distribution such as the one implied by 
Eq. (2.13). Whith these assumptions, the mean square 
relative displacement in Eq. (3.1), for an arbitrary form 
of the propagator Gij(q n ) in Eq. (2.13), is given by (here 
G(q n ) is the diagonal element Gu(q n )) ■ 



0(z)H[u 4 (z)-u 4 (O)] 2 ) 

= 2d 1 _TY,G{q n ) [l-cos(g„*)] 



(3.2) 



while the mean projected area of a given flux line is given 

by 



(u 2 ) =d 1 _Tj2G(q n 



(3.3) 



Having derived the structure factor through our ap- 
proach, we want to compare it to the corresponding quan- 
tity derived through the Boson mapping and similar ap- 
proaches. It has been shown in reference 2 using the boson 
mapping that this last quantity, in Fourier space, is given 
by: 



S(q,q z ) 



pTq 2 /n 



<Z 2 + e 2 (q)/T 2 



(3.4) 



where the excitation spectrum of the equivalent bosons 
has the usual form 



£(q) 

T 



V 2k 



2 \ 2 



pq 2 V(q) 



1/2 



(3.5) 



From equation (3.4), it is not difficult to obtain that the 
partial Fourier transform 5(q, z) has the the following 
expression 



S(d,z) 



^5(q,«,)e* 



pe 



-£(q)N/T 



^l + A P V{q)n/(Tqf 



(3.6) 



We first note that the above analytic form of the struc- 
ture factor, Eq. (3.1), and the hydrodynamic result (3.6) 
cannot be brought to agreement. In other words, no 
Gaussian approximation of the form (2.13) can produce 
a structure factor of the form (3.6). 

As we already mentioned in the introduction, it has 
already been noticed 4 ' 9 that the boson mapping gives 
nonsensical results when applied to flux liquids in the 
presence of splayed columnar disorder. In Appendix B, 
we comment on the definition of the chemical potential 
of flux lines in the boson mapping, which appears to us a 
little problematic. While this may be considered to be a 
rather "cosmetic" issue, in what follows we show that the 
boson mapping yields other unphysical results (in addi- 
tion to those aleady noticed in ref. ) even in the simplest 
case of pure, disorder free vortex liquids. To this end, 
let us make the following two observations regarding the 
analytic form of S*(q, z) in the boson scheme : 

(i) By taking the limit k — > oo in equation (3.4), one 
should be able to obtain the structure factor of a liquid 
of rigid flux lines, which we expect to be independent of 
z. Taking the above limit in equation (3.6), we see that 
the z-dependence does indeed cancel out, since e(q) — > 
as k — > oo. However, the prefactor 



S(q,z = 0) 



pTq 2 /n = p 
2(e(q)/T) y/l + 4pV(q) K /(Tqf 



also goes to zero in the limit n — > oo, which means that 
the whole structure factor S*(q, z = 0) will vanish as k 
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gets very large, which is not exactly what we would ex- 
pect for a liquid of rigid flux lines. In contrast to this 
behaviour, the structure factor derived within our vari- 
ational perturbation theory is well-behaved in the limit 
k — ► oo , and it is not difficult to see from Eqs. (3.1)-(3.3) 
that we obtain, for an arbitrary choice of Gij(q n ) (with 
either /j, = or a nonzero /i), in the above limit : 



S(q, z)=p[l + pg (r)] 



(3.7) 



which is the correct result for a liquid of rigid hard rods. 

(ii) For an "ideal gas" of non-interacting flux lines, the 
mean projected area (u 2 ) of a given flux line is given by 



(u 2 ) 



d±TL 
12k 



(3.8 



and the mean square displacement <f>(z) is given by 

2T 

4>{z) = — \z\ (3.9) 

K 

Using these expressions in Eq. (3.1) above, one finds, for 
non-interacting flux lines in the thermodynamic (L — > 
oo) limit, 



Sid.(q, z) = P exp ( - -^7 |z|) 



(3.10) 



Comparing this expression to the expression (3.6) of the 
structure factor of an interacting flux liquid in the boson 
mapping and given the expression (3.5) of the Boguli- 
ubov spectrum, we see that, at any given value of q, 
£(q, z) decays more rapidly as a function of z than its 
free counterpart Sid. fa z). This result is very surprizing, 
to say the least, as it is very difficult to imagine a flux 
line assembly where the density decorrelates on shorter 
length scales and vortices fluctuate more strongly than 
in an "ideal gas" of freely fluctuating flux lines. Indeed, 
the effect of the repulsive interactions between flux lines 
is to reduce the fluctuations of vortices, as is evidenced 
at the freezing transition: in the absence of interactions, 
a liquid of vortices would never freeze into a vortex solid 
and would continue to have large fluctuations, Eqs. (3.8)- 
(3.9), down to T = 0. Going back to equation (3.1), it 
is not difficult to verify that our result for the structure 
factor does not exhibit the above behaviour. With the 
particular choice (2.24) for the pair distribution function 
go{r), we obtain 10 



Sfaz) = p 2 3o(q)e + 

+ pexp( ^— (i-c-y^mT^N)) (3.ii) 

which has the following, nonvanishing limit as z — > oo 
Sfa z - co)=p 2 .g (q)e~37fe + pcxp ( - —2==) 

(3.12) 



indicating that the flux line densities pfa z) at both sides 
z = and z = L of a superconducting sample are not 
very much decorrelated, in agreement with the experi- 
mental findings of rcf. 22. In this reference, the authors 
have performed a two-sided decoration experiment on 
both sides of a superconducting sample, and observed 
the same pattern of vortices in each case. If the struc- 
ture factor obeyed the boson mapping result, Eq. (3.6), 
the densities /5(q, 0) and pfa L) would have been com- 
pletely decorrelated and very different from one another. 
The fact that both density patterns in the experiments 
of ref. 22 look exactly the same might be evidence of the 
fact that there is a strong correlation between pfa 0) and 
pfa, L), and that Sfa, z) is finite in the limit z — > oo. 

In fact, one can even attempt a slightly more quan- 
titative comparison of our results, Eqs. (3.11)-(3.12), 
and the experimental findings of rcf. 22. From Eqs. 
(3.11)-(3.12), it is not difficult to see that the ratio 
Sfa, z = 0)/Sfa, z = L) is of order unity for q not too 
close to 0, and is such that 



S(q,* = 0) 
Sfa,z = L) 



1 as q — > 



(3.13) 



in agreement with the decoration experiments of Yoon et 
al. It would be interesting to do a more detailed analysis 
of the experimental results in light of the present model 
(taking into account nonlocal elasticity 27 ), as such an 
analysis may lead to values of the elastic constants of the 
flux lattice (in particular, of the tilt modulus k) which 
are closer to the conventional values of these quantities 
(recall that the analysis of the experiments of ref. 22 us- 
ing the boson mapping yields elastic constants which are 
three to four orders of magnitude below what is expected 
from the standard theory of elacticity of vortices in uni- 
axial type II superconductors 27 ' 25,26 ). 



IV. CONCLUSION 

In summary, in this paper we have extended the ap- 
proach developed in (I) to better take into account the 
effect of thermal fluctuations of flux lines in vortex liq- 
uids. Using a variational Hartrec approximation, we 
have shown that thermal fluctuations strongly reduce the 
"mass" term which confines the fluctuations of the inter- 
nal modes of flux lines. We have also argued that our 
approach, which is based on the use of the conformation 
variables {r n (z)} as the true dynamical variables in terms 
of which all (Gaussian) averages are taken, yields physi- 
cally more reasonable results than the boson mapping 2-4 
or other hydrodynamic approaches 7 which, in contrast, 
use the density p(r, z) as the basic dynamical variable of 
the system. In (I) , we had argued that a careful numeri- 
cal measurement of the pair distribution function for the 
cm. mode go(r) could shed more light on the confining 
mass term ^ m a vortex liquid. Here, we have shown 
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that a numerical measurement of the temperature de- 
pendence of the mean square projected area of flux lines 
can also be used to test the predictions of the present 
approach. In particular, a nonlinear temperature depen- 
dence of (u 2 ) could be the signature of a weakly entangled 
state, where vortex fluctuations, as measured by (u 2 ), do 
not diverge with the sample thickness L. 

It is worth noting that the considerations of (I) and of 
the present work can be easily extended to the dynam- 
ics of flux liquids. This will be the subject of a future 
contribution 29 . 



V. APPENDIX A : MINIMIZATION OF THE 
VARIATIONAL FREE ENERGY 

In this Appendix, we present a few details of the min- 
imization of the free energy of Eq. (2.19) with respect 
to the propagator of internal fluctuations Gki{q m )- Tak- 
ing the derivative of Eq. (2.19) with respect to this last 
quantity amounts to computing the following expression 



dF v 



i r N 



eff 



-PH eff [{r 0i }] 



(5.1) 
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dG k i(q m ) Zi J dG i3 (q n ) 
with the effective Hamiltonian of the cm. modes: 

H eff =H -Tln<S>{g) , (5.2) 
and the partition function 

N 



Zi= f f[dr 0l e^""' 
J i=i 



[{ro,}] 



(5.3) 



The derivative of the first term in H e f f can be found in 
a straightforward way, and gives : 



dH 



dGki {q n 



= LTd ±K qf n 6 kl - \ TL V / q 2 V(<i)e^^- r ^ e"^ 2 ^ [(5 ik + 6 jk )S k i - 25 ik 5 5 i] . (5.4) 

z i.i A 



On the other hand, the derivative of the second term in Eq. (5.2) can be taken in a standard way, with the result 

d 



In $(0) = Ld ± [G~ 1 (q v 



kl 



(5.5) 



dG k i{q m ) 

With hindsight from the results of (I), we can expect (u^(z) • Uj(z)) ~ 0, so that 

dj = ([u t (z)-u J (z)] 2 )/d ± c2(u 2 )/d ± . 
Using this last approximation into Eq. (5.4), and averaging over the cm. positions in Eq. (5.1), using the fact that 



1 r N 2 f 

^ J Y[dr 0i e^ r -- r <«) e^//« r ->] = — ^ _ - J dr 0i dr 0j g (r 0t - r 0j )e^^-^ 

P 



N - 1 



.go (q) 



(5.6) 



we immediately arrive at Eq. (2.20a)-(2.20b) of the text. 

VI. APPENDIX B : STRUCTURE FACTOR OF 
THE FLUX LIQUID 



In this appendix, we show some details of the deriva- 
tion of the structure factor 



S{v,z;v',z') = {p{v,z)p{v',z')) 



(6.1) 



where the statistical average (• • •) is taken with the statis- 
tical weight exp(— f3H)/Z, with H the Hamiltonian de- 
fined in Eq. (2.11) and Z = Z Z U , with Z = Ti(c- fiiia ) 



and Z u = Tr(e~' 3ff "). We expect flux liquids at equi- 
librium to be translationally invariant, and the structure 
factor to depend only on the relative coordinates (r — r') 
and (z — z'), i.e. S(r, z; r', z') = S(r — r', z — z'). As 
a consequence, we have for the Fourier transform of the 
density-density correlation function : 

(p(q, 1z)M, q' z )) = (2^) d «5(q + q')5(q z + q' z ) S(q, q z ) 

(6.2) 

In the following, we shall be interested in the quantity 
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S(<1 <lz) — ~ 1 )p(~ Q ~<Z )) (6 3) where we used equation (6.2) above and the fact that 

LL d l ' ' 1 (27r) d (5(q = 0)<5(<? z = 0) = LL d ± in the limit L, L ± -» oo. 

We have : 

(p(q,g,)p(-q,-«»)) = | drdz | dr'dz' (p(r, z)p{v> , z')) e -i*(r-r') e - ig , {z - z >) 

N N . . 

= ^2Y1 / drdz / dr ' dz ' ( 6 ( r ~ r i( z )) 6 ( r ' ~ r i( z '))) e - 4q -( r - r ')e-^( z - z ') 
i=i j=i •> J 

= ^2Y1 I dz I dz ' (e" ,q ' (rorro,) )o(^ q ' ,u ' (z)_Ul(z ' )) }i^ ( "'>l 
»=i j^i ^ 

AT 

+ J2 dz ^'(e" iq ' (UiW " U;,(z ' )) )i e-^ (z - z ' } 

i=l ^ 

N r r 

»=i j^i ^ 

N r r 

-f ^ dz d Z 'e _ 29°W<[«.,a(^)-«.,a(z')]["»,3(^)-«i,/3(z')]>l e -^( Z ~2') (6.4) 

i=l ^ ^ 

In the above equations, (• • -)o and (• • •)„ denote averages over the cm. and internal modes, with statistical weights 
exp(— (3H ) and exp(— (3H U ) respectively (H and H u are the Hamiltonians given in equations (2.18) and (2.13) of the 
text. In our variational model, we find that the internal degrees of freedom Ui(z) and Uj(z) belonging to two different 
lines i ^ j are decoupled. We therefore can write, for i ^ j : 



([Ui,a{ Z ) ~ u iA Z ')\[ U iA Z ) ~ u 'jA z ')])l = (Ui,a{ Z ) U i,p( Z ) + u jA z ') u jA z '))l 

V> (6-5) 



d± 

where, in going from the first to the second line, we used the fact that (ui iCt UiAi — SajiAiA- Equation (6.4) becomes 

JV 



(p(q,«*)p(-q,-&)> =EE I dz I dz ' (c- l ^- ro Aoe-^ q2{u2) e^*"^ + 
i=i jjti J •> 

N f f 

+ dz dz'e~^ 

i=i •* J 



N 

q a q(}([ui, a {z)-Ui, a {z')][ui,0(z)-Ui t D(z')]) 1 p -iqz{z-z') 



c - l i^ z - z ' (6.6) 
Using the fact that J dzj dz' e ~ lqz( - z ~ z ) = L 2 5 qzfi , and noticing that 



JV 



^^ (e -iq.r oie -iq'.r OJ)o = (27r )^( q + q ') ^(q) ( 6 . 7 ) 

i=l j^i 



which gives us here (with q' = — q) 



N 

^£ (e -iqT 0ie »q.r 03 . )o = (27r )<^ (q = Q) p2 ffo (q) = L*Vso(<0 (6.8) 

i=l j^i 



we finally obtain : 



(P(q,gz)p(-q,-«z)> - £%„o Ll- p 2 . 9 o(q) e^ 9 ^ dz dz' ^^'^ (6.9) 

i=l 



where we used the fact that {[uiA z ) ~ UiA z ')][ u iA z ) — UiA z ')])u — -ff- ([ u i( z ) — Ui(z')] 2 )„. The partial Fourier 
transform of this last expression with respect to q z leads directly to expression (3.1) of the text. 
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VII. APPENDIX C : CHEMICAL POTENTIAL OF 
A FLUX LINE LIQUID 

In this appendix, we briefly comment on the "chemical 
potential" of a flux line liquid adopted in references 2,3, 8 ' 4 , 
which seems to us to be a little unconventional. In what 
follows, we shall be considering, in the notation of the 
text, a flux line liquid in a superconducting sample of 
thickness L along the z axis. We will call O the area of 
the sample in the (x, y) plane. The self energy per unit 
length of a single flux line will be denoted by s\, and is 
given by £1 = s ln(X ab /^ ab ), where X ab and £ a h are the 
London penetration depth and the coherence length in 
the basal (ab) plane. The magnitude of the magnetic in- 
duction inside the sample will be denoted by B, and it is 
related to the flux line density p — N/Cl by B = p<j)o. 

To start with, we make the observation that the chem- 
ical potential per unit length of a system of TV flux lines 
can be easily obtained using the thermodynamic identity 
(the chemical potential here should not be confused with 
the mass of the internal modes of flux lines) 



i 1 = 



1 dF 
L ON 



(7.1) 



where F is the Hclmholtz free energy of the flux line 
system. Using the fact that the free energy per unit vol- 
ume / = F I (flL) satisfies the following thermodynamic 
equality 



df_ = 

dB 

and remembering that B 
obtain 



H 

47T 



M = 



47T 



(7.2) 

n<p — N(f>o/fl, we readily 
(7.3) 



Instead of the Helmholtz free energy F, Nelson and 
Seung 2 consider the Gibbs free energy G — F — 
V(BH/Air), which for a system of flux lines is given by : 

N x , 2 

G = J2 £i / dz i 1 + ( dr */ dz ) 2 ) + 

i=l J 

+ \llj dz V ^(z) tj(z)] -V™ (7.4) 

where V(r) = 2eoK a (r/\ ab ) is the interaction potential 
between flux line elements at height z. Upon expanding 
the square root on the rhs of this last equation, and read- 
justing the tilt modulus E\ — > k = e 2 Ei to take anisotropy 
into account, (and using the fact that the volume V = flL 
and B = Nfo/Sl), they obtain : 

G = (ei — ^) NL+^njrJdz {dn/dzf + 

l—l 

+ lY,Jd*V[r i {z)-T j {z)] (7.5) 



They then identify the term proportional to N in this 
last expression, i.e. (ei — H(f> /ATr)NL with the "usual" 
chemical potential term — pN and obtain in this way the 
chemical potential per unit length as 



which does not coincide with expression (7.3). 



(7.6) 



At the origin of this disagreement is a confusion in ter- 
minology that is worth clarifying. When we use terms 
such as "chemical potential" in vortex physics, we are 
borrowing a terminology whose most natural realm is 
the theory of gases and liquids. There, one defines the 
Helmholtz free energy F by Fu q = —T\nZ, where Z is 
the canonical partition function of the system (we use 
the subscript liq to distinguish quantities pertaining to 
liquid state theory from quantities in vortex physics bear- 
ing the same name) . Fu q contains information about the 
interactions among the particles in the system and pos- 
sibly, as is the case for vortices, about their self-energies. 
One also defines a Gibbs free energy Gu q by Gn q — fj,N, 
where \i = (dFu q /dN) is the chemical potential. In the 
grand-canonical ensemble, the grand canonical potential 
Q gr = —TlnZgr can be written in terms of Fu q and Gu q 
as Q gr = F Hq - Gug = F Uq - fj,N. For homogeneous sys- 
tems, the equilibrium density of particles p is obtained 
by minimizing the quantity (Q gr /V) = f — p,p, using the 
equation 



dQ 



gr 



dp 



= 



(7.7) 



or if a one body external potential is present, using a 
functional derivative 



Sp(r) 



= 



(7.8) 



This last equality is in fact the expression of the second 
Hohenberg-Kohn-Mermin theorem 14,15 which is the basis 
of density functional theories of classical liquids. 

In the case of flux line systems, what we call Gibbs 
free energy G = F —V{BH /Air) corresponds to the grand 
canonical potential Q gr of liquid state theory, as can eas- 
ily be seen by comparing the definition of G for vortices, 
G = F-L(<p H/4w)N, to Q gr = F-pN, and remember- 
ing that the equilibrium density of flux lines p — B/<po is 
obtained in a way reminiscent of equation (7.7), i.e. 



dG_dG_ Q 
~dp~ ~dB~ 



(7.9) 



The fact that G for magnetic systems corresponds to the 
grand canonical potential has been recognized a long time 
ago by de Gennes 13 , who obtained the pressure in the flux 
line system as P = —(dG/dfl), which is reminiscent of 
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thermodynamic relation Q gc = — PV (with V the volume 
of the system) . In view of the fact that the Hclmholtz free 
energy F has the same meaning in liquid state theory and 
for magnetic systems, one is led to identify V(BH /At:) 
with the Gibbs free energy fiN of liquid state theory, and 
this identification leads to the result (7.3) for the chemi- 
cal potential per unit length of the flux line system (see 
also the table below, which contains a summary of the 
correspondence between flux line and liquid quantities) 

We now can understand the problem with the deriva- 
tion of the chemical potential \x in reference 2 . The prob- 
lem with this derivation is that the term e\NL belongs 
to the Hclmholtz free energy F while the true chemical 
potential should be identified only form the Gu q part of 
the grand canonical potential, i.e. from V(BH/4ir) — 
L(<fioH/4:ir)N. This improper identification of the chem- 
ical potential p is, as far as the results of refs. 2-4 are 
concerned, only a minor, "cosmetic" point. However, in 
order to obtain meaningful results in a grand-canonical 
formulation (using an approach similar to ours in (I) and 
in the present work) , it is very important that the correct 
expression of pi, equation (7.3), is properly identified. 



Flux lines 


Classical Fluids 


Helmholtz free 


Hclmholtz free 


energy 


energy 


F = -TlnZ 


F = -TlnZ 


Gibbs free energy 


Grand canonical 




potential 


G = F — L(<l> H/4ir)N 


Q gr = F - pN 


(dG/dB) = 


(dQ g r/dp) = 


chemical potential 


chemical potential 


per unit length 




ix = (df/d P )/L 


M = (df/dp) 







Table 1. Summary of correspondence between vortex 
and classical fluids quantities. 
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